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Small-scale magnetic fields on late-type M-dwarfs 
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Abstract. We performed kinematic studies of the evolution of small-scale magnetic fields in the surface layers of M-dwarfs. 
We solved the induction equation for a prescribed velocity field, magnetic Reynolds number Rbm, and boundary conditions 
in a Cartesian box, representing a volume comprising the optically thin stellar atmosphere and the uppermost part of the 
optically thick convective envelope. The velocity field is spatially and temporally variable, and stems from detailed radiation- 
hydrodynamics simulations of convective flows in a proto-typical late-type M-dwarf (T e g = 2800 K, log g — 5.0, solar 
chemical composition, spectral type ~M6). We find dynamo action for ReM > 400. Growth time scales of the magnetic 
field are comparable to the convective turn-over time scale (« 150 sec). The convective velocity field concentrates the 
magnetic field in sheets and tubular structures in the inter-granular downflows. Scaling from solar conditions suggests that 
field strengths as high as 20 kG might be reached locally. Perhaps surprisingly, ReM is of order unity in the surface layers 
of cooler M-dwarfs, rendering the dynamo inoperative. In all studied cases we find a rather low spatial filling factor of the 
magnetic field. 
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^ ■ 1. Introduction 



M-type dwarfs show the highest degree of magnetic activity. 
Besides X-ray and H a emission as indirect tracers, Zeeman 
broadening of magnetically sensitive photospheric lines has 
been observed (Johns-Krull & Valenti 1996, Kochukhov et al. 
2001). The photospheric lines show no rotational modulation 
or net polarization, indicating that fields of small scale (rel- 
ative to the star's radius) exist on the surfaces of M-dwarfs. 
These fields can have an appreciable strength (fewkG), and 
can cover a substantial area (filling factors up to 50 %). 

Little is known about the structure of magnetic fields in 
the photospheres of M-dwarfs. From an observational point 
of view one would like to get some input from theory, that 
would alleviate the problem of disentangling field strength 
and filling factor. To this end we started to investigate the 
kinematic effect of the convective velocity field on a magnetic 
seed field. This has become possible due to recent progress in 
the hydrodynamical modeling of atmospheres in the regime 
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of cool M-dwarfs (Ludwig, Allard, & Hauschildt 2002). The 
present investigation is targeting at an application to spec- 
troscopy, meaning that we focus on layers near optical depth 
unity. To overcome partially the limits stemming from the 
restriction to the kinematic magneto-hydrodynamics regime 
(hereafter "kinematic MHD"), we tried to relate our results 
to the situation in the Sun, emphasizing the scaling from fa- 
miliar solar conditions to M-dwarfs. 



2. Model setup 

2.1. The convective velocity 

The velocity field was taken from a radiation-hydrodynamics 
simulation of a prototypical M-dwarf atmosphere (Ludwig et 
a. 2002, see Figg) of T cS = 2800 K, \ogg = 5.0, and so- 
lar chemical composition (corresponding to a spectral type of 
wM6). The 3D radiation-hydrodynamics code of Nordlund 
& Stein (Stein & Nordlund 1998, and references therein) has 
been adapted for this purpose. The code treats the interac- 
tion of radiation and gas flows in detail, rotation and mag- 
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Fig. 1. Snapshot of the vertical component of the magnetic 
field at optical depth unity (contours), emergent intensity 
(grey tones), and the velocity field (arrows) towards the end 
of an evolutionary sequence of an unquenched model with 
Re M =20. 



netic fields are neglected. Rotation has only a minor effect on 
the structure of the convective flows near the stellar surface 
due to the short convective turn-over time scales (relative to 
typical rotation periods of several hours or longer). Magnetic 
fields can potentially influence the flow structure significantly 
provided the field is strong, i.e. dynamically important. How- 
ever, this corresponds to the full MHD case which is beyond 
the scope of this paper. We selected a temporal sequence to- 
wards the end of the simulation run, which comprises 1500 s 
of stellar time, corresponding to roughly 10 convective turn- 
over time scales. The sequence consists of 150 snapshots of 
the flow field, each comprising 125 x 125 x 82 grid points 
corresponding to 250 x 250 x 87km 3 in geometrical size. 
At any instant in time about 10 granular cells were present 
in the computational domain, ensuring a statistically repre- 
sentative ensemble. M-dwarf granulation is on the qualitative 
level similar to the solar granulation. Quantitatively, there are 
differences which are discussed below where relevant. 

2.2. Kinematic magneto-hydrodynamics 

We assume the kinematic regime of MHD, where one ne- 
glects the back-reaction of the magnetic field on the fluid 
motions. This means that the non-linear behavior cannot be 
rendered correctly in this regime, which is also sometimes 
referred to as the linear regime of MHD. Then solving the 
MHD equations reduces to the problem of seeking the solu- 
tion to the time-dependent induction equation: 

SB 

~dt 



V x (u x B) +?7V 2 B, 



(1) 



magnetic diffusivity. We assume a spatially constant mag- 
netic diffusivity. Note, that the problem is linear on this level 
of approximation. We vary the magnetic diffusion by set- 
ting the magnetic Reynolds number ReM which is defined 
as ReM = U£/i], where U is a characteristic velocity, £ 
a characteristic length scale, and 77 the magnetic diffusiv- 
ity, which in turn is related to the electrical conductivity a 
through X] = 1/ /zcr. 

We solve Eq. (]lj) using staggered variables on the grid 
of the hydrodynamical flow field. The numerical method was 
originally developed by Galsgaard and others (Galsgaard & 
Nordlund 1997 1 for general MHD purposes. A special ver- 



sion is the code used by Archontis & Dorch (Dorch 2000 



Archontis, Dorch, & Nordlund 2002) to study dynamo action 
in prescribed flows. 

We kept the same velocity field, and studied different 
magnetic diffusivities. To first order we hope to represent the 
different conditions encountered in various M-dwarf atmo- 
spheres by the the change of magnetic diffusivities. 

2.3. Magnetic quenching of the velocity field 

The kinematic approximation to the MHD equations — i.e. 
neglecting the back-reaction of the magnetic field by omit- 
ting the Lorentz force in the equation of motion — be- 
comes inaccurate when the magnetic field becomes suffi- 
ciently strong. That is, when the magnetic energy density ap- 
proaches equipartition with the kinetic energy density, cor- 
responding to a field strength B oq = u^JJTp. One can try to 
capture the effect of a near-equipartition magnetic field, by 
utilizing a quenched velocity field that has a reduced ampli- 
tude in the regions where the magnetic field is strong. In a 
few cases we employed such a quenched velocity field when 
solving the induction equation Eq. (|lj). The quenched velocity 
field u q was modeled according to 



u q = u cxp-a(eM/eK) 2 



(2) 



where u is the prescribed velocity field (from the radiation- 
hydrodynamic simulation), B the magnetic field, and r\ the 



where a is a constant, e« the local magnetic energy density, 
ex the average kinetic energy density, taken to be constant 
in time and space. The quenching reduces the velocity am- 
plitude at locations where the magnetic energy (or field) be- 
comes comparable to the kinetic energy of the fluid motions: 
it thus reduces the growth of the magnetic field in these re- 
gions causing it to saturate. We chose a = 1.75 which en- 
sures u q « u up to field strength of 10% of B cq , and u q w 
at a field strength of 1.2 times equipartition. Note, that this 
quenching procedure makes the problem non-linear in B, and 
that a constant ex merely introduces a scaling factor in the 
solution, i.e. its concrete value is unimportant. In future work 
we plan to replace Eq. (Q) with a more realistic quenching 
expression, taking into account the geometry of the flow and 
the magnetic field. 

2.4. Initial and boundary conditions 

The numerical method used to solve Eq. ([!]) by default 
considers periodic boundary conditions in all three dimen- 
sions. We implement non-periodic vertical boundary condi- 
tions through "ghost zones" at the top and bottom of the do- 
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main: then one may constrain the magnetic field to be ver- 
tical (horizontal) by requiring its horizontal (vertical) com- 
ponents) to be anti-symmetric across the boundary, and the 
remaining component(s) to be symmetric. 

The initial magnetic field is either uni-directional (vertical 
or horizontal) or smoothly varying on small scales. The initial 
field strength B is set to be much smaller than B cq (only 
relevant for the cases with quenched velocity field). 

3. Results 

3.1. Magnetic Reynolds numbers of order unity in cooler 
M-dwarfs 

In the following we discuss kinematic MHD models with 
ReM as low as 20. In the astrophysical context, low values 
of Rcm are uncommon, mostly due to the large spatial scales 
usually involved. However, in M-dwarf atmospheres we can 
be confronted with the situation of rather low Rcm (Meyer 
& Meyer-Hofmeister 1999): Fig. ^ shows the ratio of elec- 
tron to gas pressure, and ReM in the T e g = 2800 K M-dwarf 
model as well as the Sun. ReM has been evaluated assum- 
ing constant typical length and velocity scales independent of 
deptff]. ReM primarily reflects the run of the electric conduc- 
tivity in the atmosphere, which in turn is mostly controlled by 
the electron to gas pressure. The conductivity has been eval- 
uated assuming a weakly ionized plasma (Kopecky 1957, as 
quoted by Stix 1989). 

Even if one concedes an appreciable degree of uncertainty 
related to the choice of scales, at sufficiently cool tempera- 
tures Rcm reaches order unity in the surface layers. This is a 
consequence of the declining electron density, the shrinking 
of spatial scales (due to smaller pressure scale heights), and 
smaller convective velocities (due to lower energy fluxes and 
higher atmospheric densities). This refers to the surface lay- 
ers. Qualitatively, we expect a strong increase of ReM with 
depth, and beyond a certain depth the regime of ReM > 1 
is reached again. However, gas motions in this depth will 
generally be slower and the tangling of magnetic field lines 
less rapid, which may reduce the efficiency of chromospheric 
and coronal heating. Whether this plays a role for the ob- 
served decline of stellar activity at the transition from M- to 
L-dwarfs (Gizis et al. 2000) is presently a matter of debate 
(Mohanty et al. 2002, Berger 2002). 

3.2. Scaling of the field strength 

In order to obtain an estimate of the field strength to be ex- 
pected in M-dwarf atmospheres — provided local conditions 
indeed control its value — we compared kinetic energy den- 
sities and pressures obtained from our M-dwarf radiation- 
hydrodynamics model with a similar model computed for 
the Sun. Figure || shows that the kinetic energy densities in 
the surface layer of our reference M-dwarf model are sys- 
tematically smaller than in the Sun, but only within a factor 
of 2. At the t = 1 level the pressure is higher by a factor 

1 Sun: L = 1500 km, v = 2.4km/s; Mdwarf: £ = 80 km, v = 
0.16km/s 



600 




IOg 10 Xr oss 

Fig. 3. Kinetic energy densities expressed in terms of the cor- 
responding equipartition field strength as a function of Rosse- 
land optical depth. Depicted are results from hydrodynamical 
simulations for the T cS = 2800 K M-dwarf (solid lines) and 
the Sun (dash-dotted lines). The total kinetic energy density 
is broken into contribution of horizontal (thick lines) and ver- 
tical motions (thin lines). 

of 12 and density by a factor of 42 in the M-dwarf model 
relative to a solar model (not shown). Assuming a scaling 
of the magnetic field with gas pressure, we expect a maxi- 
mum field strength in flux concentrations of about 10 times 
the solar value, or about 20 kG if one scales values found in 
magneto-hydrodynamical simulations for the Sun (Stein et al. 
1999). However, the kinetic energy densities point towards a 
field of slightly sub-solar field strengths. 

3.3. Dynamo action and magnetic Reynolds number 

A flow is a fast (kinematic) dynamo when the exponential 
growth rate 7 is positive, where 

7 = lim 7( , = lim lim log(E M (i)/E M (0))/t, (3) 

?;^0 77— >0 t—>oo 

where Em is the total magnetic energy, and j v is the growth 
rate at a specific diffusivity, corresponding to a certain ReM- 
That is, fast dynamo action requires a continuous increase 
of magnetic energy, even in the limit of vanishing diffusiv- 
ity. This limit is relevant because most astrophysical sys- 
tems have Rcm 3> 1 and small but non-zero 77. It is be- 
lieved (but not proven) that turbulent astrophysical systems 
are fast dynamos operating at ReM ^> 1. When ReM in- 
creases, the length scale of magnetic islands decreases and 
_ 1 

scales as Rc M 2 . There is a maximum value of Rcm, which 
can be achieved with our numerical resolution Air = 2 km, 
and therefore the largest magnetic Reynolds number that we 
can allow is of the order of 400 corresponding to 2 Aa; (the 
Nyquist wavelength). Figure ^ shows the two-dimensional 
power spectrum of the total magnetic field and its vertical 
component at optical depth unity in the case ReM = 100: the 
power on small-scales close to the resolution limit is much 
less than on the granulation scale. Hence in the following, we 
concentrate on situations with ReM < 400. 
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Fig. 2. The ratio of electron to gas pressure (left panel), and the magnetic Reynolds number (right panel) as a function of 
Rosseland optical depth in a T e g = 2800 K M-dwarf (solid line) and a solar model atmosphere (dashed line). Note the low 
magnetic Reynolds number in the M-dwarf model, primarily reflecting the low electric conductivity of the stellar gas in the 
rather cool M-dwarf atmosphere. 
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Fig. 4. The power of the magnetic field at the surface (r = 1) 
in a model with ReM= 100. Shown are the total field strength 
(full curve) and the vertical field component (dashed). The 
vertical lines indicate the Nyquist wavelength (left) and the 
typical size of the granular cells (right). 
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Fig. 5. The total magnetic energy Em/Eo as a function of 
time in units of 150 seconds (the typical convective turn-over 
time). Five models are presented, with different amounts of 
diffusion: Re M = 800 (full), 400 (dotted), 200 (dashed), 100 
(dashed dotted), and 20 (dashed triple dotted curve). 



3.4. Kinematic models 

In the absence of non-linear effects, in case of a dynamo, 
one expects a continuing exponential growth of Em - Figure |] 
shows the result in terms of Em(£) for five different mod- 
els corresponding to varying Kgm- most of the models in fact 
are dominated by decaying modes, with negative growth rates 
•jrj < 0. In particular, the case with Rcm= 20 (high diffusion) 
is clearly an example of an anti-dynamo: the diffusion works 
faster than the flows can sweep up the field and concentrate it 
in the inter-granular downdraft lanes, and the dominant mag- 
netic mode is a decaying one. Increasing ReM decreases j v 
(numerically), so that for Rgm= 200 the decay time is 40 time 
longer than at ReM= 20. In terms of providing dynamo ac- 
tion, the most promising cases are those with lower diffusion 
and ReM> 300-400: at ReM= 400 a growing mode seems to 
be dominating. One of the models in Fig. |^ had Rcm= 800, 
i.e. it produced magnetic structures too small to be resolved, 



and therefore magnetic energy is lost from the computational 
domain and Em seems to saturate because of enhanced nu- 
merical diffusion. 

In the high diffusion case (e.g. ReM= 20) the magnetic 
field varies smoothly across the domain, and its power peaks 
on scales approximately I ps 50 km when considering B z (see 
Fig |l[). At higher ReM more small scales are generated mostly 
around the downdraft lanes. 

3.5. Quenched models 

To examine the structure of the magnetic field in the high 
ReM cases with dynamo action, we computed a few models 
including the quenching mentioned above: quenching the ve- 
locity field only affects those models that are dominated by 
growing modes (7, ; > 0). In case of a low ReM, the results 
remain unchanged since the field never becomes compara- 
ble to B oq . For Rcm= 200 the dominant mode is still decay- 
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Fig. 6. The total magnetic energy Em/E as a function of 
time for a model with ReM= 400 and a quenched velocity 
field (full). Also shown is the total available kinetic energy 
(dashed curve). 



ing. At ReM= 400 the field is amplified until it approaches 
equipartition and the quenching causes Em to saturate; see 
Fig. |[ At the time of saturation there is a lot of power at 
small scales close to the resolution limit, but the main power 
is concentrated at larger scales around 100 km. The saturation 
level is artificially fixed by our choice of quenching function 
Eq. (||), the value of a, and our choice of a constant value of 
eK- Hence our simulation results in no direct information on 
the amplitude of the magnetic field, only on the geometry that 
one may expect in a non-linear simulation. 

Figure [7] shows the PDFs of B z for three different values 
of Rejvi: the broad distribution of field strengths for low ReM 
should be seen in contrast to the narrow one for high ReM- 
At any time the strongest part of the field occupies only a 
very small part of the volume no matter what the level of 
diffusion is, but for low ReM the preference for strong fields 
is less pronounced; the field is less intermittent. We attribute 
the non-monotonic behavior of the models displayed in Fig. ^ 
to the model with ReM = 400 showing dynamo action while 
the others do not. 



4. Conclusions 

- The geometry of small-scale (i.e. locally generated or 
modulated) magnetic fields looks similar to the situation 
for the Sun (i.e. mostly located in the intergranular lanes). 
The basic reason is that M-dwarf granulation is qualita- 
tively similar to solar granulation. 

- There are differences due to potentially quite different 
magnetic diffusivities (and values of ReM)- 

- Depending on Rcm we get or do not get local dynamo 
action. 

- Since we are making a kinematic study, the absolute field 
strength is undefined. Scaling from the Sun would suggest 
flux tubes with up to factor 10 higher fields (assuming 
equipartition with the gas pressure), or slightly sub-solar 
(within a factor 2) fields (assuming equipartition with the 
kinetic energy density). 
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Fig. 7. The probability distribution function (PDF) of the ver- 
tical field component around optical depth unity. Lines show 
the distribution in unquenched models for magnetic Reynolds 
numbers ReM = 20, 400, 200 "outside-in" of the hat shaped 
PDFs. Note the non-monotonic behaviour with ReM- 
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